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We have used exact diagonalization and quantum Monte Carlo methods to study the one- 
dimensional t-J model including the three-site hopping term derived from the strong coupling 
limit of the Hubbard model. The three-site term may be important to superconducting correlations 
since it allows direct hopping of local singlet electron pairs. The phase diagram is determined for 
several values of the strength of the three-site term and compared with that of the t-J and Hubbard 
models. Phase separation, which exists in the t-J model is suppressed. In the low electron density 
region the formation of local singlet electron pairs is enhanced, leading to stronger superconducting 
correlations even for values J ft < 2. A large spin gap region extends from low electron densities up 
to high densities. In the low hole density region the superconducting correlations are suppressed at 
J ft > 2.8 in spite of enhanced pair formation. This is because the three-site term, while enhancing 

the formation of electron pairs, leads to a repulsion between holes. 

74.20.Mn, 71.27.+a, 75.10.Jm, 74.25.Dw 
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I. INTRODUCTION 

Since the discovery of high-T c superconductors, 1 the 
investigation of strongly correlated fermion systems is 
a central theme of theoretical solid state physics. It 
has been established that the relevant physics of these 
cuprate superconductors is contained in the two dimen- 
sional (2D) CuC>2 planes. Anderson suggested that 
a single-band Hubbard model would be the minimal 
model needed to describe the low energy properties. 2 The 
Hamiltonian can be written as 



H, 



Hubbard 



E 



(1) 



where denotes a pair of nearest neighbors, c\ a and 

Cj CT are the usual electron creation and destruction opera- 
tor at the site i with spin a and ri{ a = c\ a Ci a . The on-site 
Coulomb term U describes strongly correlated electrons 
in the limit U/t — > oo as well as weakly coupled electrons 
for U/t < 1. 

Strong correlations make the limit U/t — > oo especially 
interesting but also difficult to investigate. The system 
avoids doubly occupied sites in this limit, since they cost 
a large energy U . In addition to the local constraint 
of excluding double occupancy, interactions take place 
through virtual double occupancy. By a canonical trans- 
formation we get an effective Hamiltonian 3,4 which reads 
to second order in t/U : 



J-ts 



ia c ja +h.c. 



(2) 
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This model which we call t-J-ts model is the subject of 
this paper. We set J = At 2 /U and t 3 = aJ/4. The 
value a = 1 is obtained by the canonical transformation, 
but we will also consider smaller values of a to provide 
a continuous switch on of the pair hopping terms. The 
brackets (i, j, k) denote a three site term, with i ^ k being 
neighbors of j. The pseudo-electron creation operator 
c f =c 
site, ni 



_ CT ) allows only single occupancy of each 

"if + n il an d Si denotes the spin operator. 
Usually the pair hopping terms are neglected, 4,5 and we 
obtain the familiar t-J model for a = 0. We note that it 
has been shown by Zhang and Rice that the t-J model 
can also be derived from a three-band model, 6 describing 
the CuC>2 planes. 

We can write the interaction term in H t -j-t 3 alterna- 
tively as a number operator and a hopping term for sin- 
glet pairs of electrons: 



± a j]T(i^ + i + h.c.), (3) 



where P) 



(1/V2)(c: t c: +U 



c ll c ]+i-f) is the creation 
operator of a nearest neighbor singlet pair. This allows 
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electron pairs to hop without the cost of breaking a bond 
J. We expect therefore at least in the low electron density 
region a strong stabilization of local electron singlet pairs. 
If this stabilization is strong enough, we will get a short 
range resonating valence bond (RVB) state, forming a 
gas of local singlet pairs. The lowest spin excitation of 
this short range RVB state will then be the excitation 
of a local singlet pair to a triplet state, and hence have 
a finite spin gap. An important difference to the usual 
t-J model (a = 0) is that the local singlet electron pairs 
of this short range RVB phase can gain kinetic energy 
of the order of a J due to the pair hopping term and we 
may expect an expansion of this spin gap phase to higher 
electron densities. Since a finite spin gap is characteristic 
of this state, we will also investigate the region of a finite 
spin gap the t-J-ts model in this paper. 

Also related to the kinetic energy is the occurrence of 
phase separation. The actual structure of the electron 
rich and poor phases is determined by a compromise of 
a gain of kinetic energy and a cost of exchange energy, 
and therefore depends on the ratio J ft and the lattice 
structure. In the one-dimensional (ID) t-J model, the 
phase separation starts at J ft ~ 2.8 (for small electron 
density) and 3.6 (near half filling). 7,8 On the other hand, 
the critical J ft value for phase separation is smaller near 
half filling in two dimensions. 9 The phase separation in 
the t-J-ts model would be quite different from the t-J 
model. This is because even for large J ft the system now 
can gain a large kinetic energy by pair hopping and con- 
sequently the phase separation is expected to be strongly 
suppressed. 

To study the effects of this new pair hopping term in 
detail, we have limited our calculations to the ID case 
for simplicity. 



II. NUMERICAL METHODS 

Two complementary numerical algorithms have been 
used to investigate the properties of these systems, ex- 
act diagonalization of small clusters by the Lanczos 
algorithm 10,11 and the Quantum Monte Carlo (QMC) 
world line algorithm. 12 

A careful choice of the boundary conditions is impor- 
tant for exact diagonalization. One reasonable choice is 
the boundary condition where all one electron orbitals 
in the noninteracting case are either fully occupied or 
empty (closed shell boundary conditions, CSBC). 7,13,14 
For ID single-band models, this corresponds to periodic 
boundary conditions (PBC) for systems with N = Am + 2 
electrons and antiperiodic boundary conditions (ABPC) 
for N = Am, with m being an integer. With this bound- 
ary condition the ground state is a spin singlet. Alterna- 
tively choosing the opposite boundary conditions, APBC 
for N = Am + 2 and PBC for N = Am will be called 
open shell boundary conditions (OSBC). If not otherwise 
stated, CSBC is used for all calculations. 



For the QMC calculations, a parallel version 15 of 
the world line algorithm with a four site cluster 
decomposition 16 has been used. The advantage of QMC 
methods is that they allow simulations of considerably 
larger systems. We have calculated systems of up to 64 
sites at an inverse temperature of fit = 32, since the 
negative sign problem did not appear for these models. 
At such low temperatures, the measured quantities do 
not change qualitatively compared to zero temperature. 
The simulations have been restricted to the subspace of 
zero winding number. 8 To control the systematic error 
of order 0(Ar 2 ) due to the finite Trotter time step At, 
measurements have been performed at Art = 0.25 and 
0.1 and extrapolated to Art = 0. 



III. PHASE SEPARATION 

The ground state of strongly correlated electron sys- 
tems is not always homogeneous, but sometimes is phase 
separated, consisting of a mixture of two phases. The 
phase separation line divides the region of stable homo- 
geneous phases in the parameter space from the mixed 
phases region. Whereas the repulsive Hubbard model is 
homogeneous for the whole parameter range, 17-19 the t-J 
model is phase separated for large J/t. 7 ' 8 

The onset of phase separation is estimated as a di- 
vergence of the compressibility. We used a finite size 
approximation 



L 



# 2 E (L; N + 2) + E (L; N - 2) - 2E (L; N) ' 



(4) 



with Eo(L; N) being the ground state energy of a system 
with N electrons in L sites. 

For the low electron density region, the phase sepa- 
ration boundary can be estimated by an argument of 
Emery, Kivelson, and Lin. 20 We start with a Heisenberg 
chain of length # > 1 in a phase separated state. Now 
we assume that two electrons evaporate from the chain 
to form a singlet bound state at a large distance from 
the chain, and compare the energy of this state with the 
initial state energy. Solving the two electron problem, it 
is found that a singlet bound state appears if 



J > J D 



2t 



1 



and the energy of this state is 



At 2 



palr = -(1 + a)J - ^—^ 



(5) 
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Evaporating the two electrons costs the energy of two 
bonds, 2Jln2. Phase separation starts when impair > 
— 2 J In 2 or at J > J cr with 



Jc 
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FIG. 1. Boundary line of phase separation in the t-J-ts 
model for different values of a, calculated by exact diagonal- 
ization of a lattice of length L = 16 and CSBC. The phase 
separation line is shifted towards larger J/t with increasing 
the pair hopping term and disappears completely. 



in the low electron density limit. This argument as- 
sumes J cr > Jp a i r , which is always fulfilled. The esti- 
mate for J cr would be exact if there are no bound states 
with more than two electrons. Following this argument, 
phase separation disappears at low electron densities for 
a > a CT = 2 In 2 — 1 0.386. In Fig. 1 we have plotted 
the phase separation boundaries in the phase diagram 
defined in the plane of the electron density n = N / L 
and the coupling J/t for different values of a. The fig- 
ure shows that at a = 0.4 the phase separation region 
survives only near n = 1 in agreement with the above 
argument. We show only two points for the phase sep- 
aration line at a = 0.4, since here the phase separation 
line diverges very rapidly to infinity for n < 3/4. 

We can see in Fig. 1 that also in the higher electron 
density region the onset of phase separation is pushed 
away towards larger values of J/t with increasing a. At 
a = 1 phase separation completely disappears for all den- 
sities n. We can see in Fig. 4(b) that the compressibility 
k is strongly reduced for a = 1 in contrast to the diver- 
gence in k for a = and 0.4. This point will be discussed 
in more detail in Sec. V. Very small values of the com- 
pressibility k in the low hole density region may indicate a 
repulsion among holes. We can see this more explicitly by 
by considering the probability of finding a single hole be- 
tween two electrons, p s ^ = (n 8 '/}j'_|_in 8 '_|_2}/(l — n). A value 
of p s h close to one indicates that practically all holes are 
isolated, while p s ^ — > suggests that the holes attract 
each other. The results in the low hole density region 
shown in Table I indeed show that p s ^ becomes close to 
one for a = 1, and demonstrate a repulsion among holes. 
This effective repulsion is produced by the large gain of 
energy of the order of aj due to pair hopping, rather 



TABLE I. Three site correlation of two particles and one 
hole p s h = {n l h l j r \n l j r 2) / (1 — n) for a density n = 7/8 in the 
t-J-ts model, calculated by exact diagonalization on a lattice 
of length L = I 6. 
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0.962 


0.947 


0.923 



than the strongly renormalized small direct hopping. In 
the low hole density region, pair hopping enhances prop- 
agation of single holes, because it can gain more kinetic 
energy than propagation of bound hole pairs, and hence 
produces an effective repulsion among holes. This strong 
repulsion is the reason of suppression of phase separation. 

Attention has to be paid on the shape of the phase 
separation line. The Maxwell construction shows that 
the phase separated region consists of a mixture of the 
two homogeneous phases that are on the crossing points 
of the phase separation line and the line of constant cou- 
pling J/t. In t-J model, the phase separation line crosses 
the line of constant J/t only once. In this case, the elec- 
tron poor phase always is a vacuum. In the t-J-ts model, 
the phase separation line crosses the line of constant J/t 
twice at certain values of a and J/t. Then the electron 
poor part is a dilute gas of bound pairs and the electron 
rich phase is a Heisenberg chain with a finite density of 
holes inside. 



IV. SPIN GAP 

The spin gap is the excitation energy from the 
singlet ground state to the lowest triplet state 
A s = lmu-yoo A S (L; N = Ln), where A S (L;N) = 
E (L; N; S z = 1) - E (L;N;S Z = 0). The determina- 
tion of the spin-gap region is an important issue in the 
characterization of a model with spin degrees of freedom. 
As it was stated already in the introduction, the presence 
of a spin gap can indicate a short-range RVB mechanism. 

A finite spin gap is present in several extended t-J 
models. One example is the t-J- J' model, where a next 
nearest neighbor Heisenberg interaction J' is added. It 
has a finite spin gap in the low hole doping region. 21 Elec- 
tron pairing is also favored by adding a repulsive next 
nearest neighbor interaction V and this effect introduces 
a finite spin gap in the t-J-V model at quarter filling. 14 
One common characteristic of these models is the restric- 
tion of the spin gap to the vicinity of the phase separation 
line. A small spin gap region exists also for the t-J model 
in the low electron density region at 2.0 < J/t < 2. 95. 7 

The spin gap region for the t-J-ts model is determined 
for densities of n = 1/2, 2/3 and 4/5 by exact diagonal- 
ization of chains with different sizes. One needs to use 
CSBC and OSBC as explained in Sec. II. We have used 
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FIG. 2. Extrapolated value of the spin gap for the t-J -is, 
model with a = 1 at several densities as a function of the 
coupling J/t. 



the same fitting function as Refs. 14 and 21 for extrapo- 
lation to L — > oo: 



A? C (L) = A s 



,BC 

h 

L 



,BC 



_L 2 



where BC denotes the boundary condition. 

In the low electron density region, the spin gap region 
is estimated to start where a singlet bound state of elec- 
tron pairs is formed, that is at J = J pa ; r (Eq. (5)). With 
increasing a, J pa ; r decreases, and the onset of the spin 
gap region shifts therefore towards smaller couplings J/t 
for n <C 1- For a = 0.4 and 1 we have found numeri- 
cally that the spin gap region extends up to densities of 
n > 2/3. Also for these densities, we have found that 
the spin gap region starts at smaller couplings J/t with 
increasing a. The value of the spin gap for a = 1 can be 
seen in Fig. 2 for different densities n as a function of the 
coupling J/t. 



V. CRITICAL EXPONENTS 

In ID, we can describe the long range fluctua- 
tions of systems belonging to either of universality 
classes Tomonaga-Luttinger liquids (TLL) or Luther- 
Emery liquids (LEL) 22 by a single correlation expo- 
nent K p . 18 i 19 i 22 ~ 25 Tomonaga-Luttinger liquids are char- 
acterized by gapless charge and spin excitations. The 
long range correlation functions of TLL decay as power 
laws. The ID Hubbard model falls into this universality 



eh 



and K p decreases for all densities n ^ 1 from 



K p = 1 in the free electron limit U/t — > to K p = 0.5 
in the strong coupling limit U/t — > oo. Another example 
is the t-J model before phase separation. 7,8,25 Models 



with a finite spin gap and have gapless charge excita- 
tions generally and belong to the LEL universality class. 
Charge density wave (CDW) and singlet superconducting 
correlations have power law behavior, spin density wave 
(SDW) and triplet superconducting correlations decay 
exponentially because of the spin gap. Superconduct- 
ing correlations dominate the long range fluctuations for 



both universality classes if K p > 1. This 



for the 



t-J model in a narrow region of the phase diagram before 
phase separation. 7,8 

It is known from conformal field theory that K p can be 
calculated from numerical results of the low energy spec- 
trum in finite systems alone. The results thus obtained 
have to be checked for consistency a posteriori. We cal- 
culate K p from a combination of two relations. The first 
one expresses K p by the compressibility k and the charge 
velocity v c : 23 



K n 



(9) 



Since the charge velocity v c is difficult to determine from 
the spectrum, we have used a second relation to eliminate 
v c . The Drude weight <7o of the ac conductivity is given 
by the energy shift of the ground state in the presence 
of a field 26 and for ID systems also by the correlation 



(8) exponent 



Co 



ttL 



2v,K 



(10) 



</>=o 



where Eo(<f>) denotes the ground state energy of the sys- 
tem with twisted boundary conditions with a phase factor 
<f>. Together with Eq. (9) we get 



K, 



-aon 2 K. 



(11) 



The phase diagrams for a = 0.4 and 1 are shown 
Fig. 3, and we discuss them by comparing with the phase 
diagram for the t-J in Ref. 7. In the limit J/t — > the 
correlation exponent is K p = 1/2 for any a and n ^ 1 as 
in the same limit of the t-J model. Before phase sepa- 
ration, the value of K p increases monotonically with the 
coupling J/t, and we find as a precursor to phase sep- 
aration dominant superconducting correlations K p > 1. 
The reason for the strong increase of K p before phase 
separation is the divergence in the compressibility k , as 
can be seen in Fig. 4 where we show along with K p the 
compressibility k and the Drude weight oq. Since there 
is no phase separation for a = 1, the superconducting 
region extends up to J/t — > oo and the divergence in the 
compressibility k and in K p is also removed. 

The effect of the pair hopping term is different depend- 
ing on the electron density n and the coupling J/t. In 
the low electron density region we have a strong enhance- 
ment of local singlet electron pairs, which can move by 
the pair hopping term without the cost of breaking a 
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FIG. 3. Phase diagram of the t-J-ts model from L = 16 
cluster calculation, (a) a = 0.4, and (b) a = 1. The solid 
curves show the contours of constant correlation exponent K p , 
calculated from for systems with an even number of electrons 
and CSBC. The dashed-dotted line represents the border line 
of phase separation, the dashed line the border line of the 
region of finite spin gap. 



bond J. These mobile local electron pairs enhance su- 
perconducting correlations and K p is increased with a 
for all densities n at low couplings J/t < 1. On the 
other hand, we have a strong suppression of K p with in- 
creasing a for low hole densities and couplings J/t > 1 
and a slight suppression of K p for low electron densities 
and large J/t. The reason for this suppression of K p 
is that the enhanced propagation of single holes in the 
low hole density region discussed in Sec. III. This en- 
hancement produces an effective repulsion among holes, 
which reduces drastically the compressibility k, as seen 
in Fig. 4(b). Therefore the value of the correlation ex- 
ponent K p via the relation Eq.(ll) is also reduced, even 
though the Drude weight <7o increases with larger values 
of a. 

In addition, we have also shown the border line of the 
spin gap region in the phase diagrams Fig. 3. The contour 
lines of the correlation exponent seem to be continuous 
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FIG. 4. (a) Correlation exponent K p , (b) compressibility 
k, and (c) Drude weight ao for a = 0, 0.4 and 1 as a function 
of J/t at n = 3/4. Calculated by exact diagonalization of a 
lattice of length L = 16 and CSBC. 



at this transitions line, although the long range behavior 
of the correlation functions changes drastically. 



VI. CHARGE AND SPIN STRUCTURES 

We examine the charge and spin structure factors in 
more details and calculate them by a QMC method. The 
charge and spin structure factors are defined as: 

1 L 

5'charge(fc) = - ^ e lk(j m ) {{n jA + n jtX ) (n m>t + n m ,j.)) , 

j,m 
1 L 

S spm (k) = -J2e lk(j m) ((nj,t - nj,l)(nm,t ~ n m ,i)), 



L 



_£y*«- m )<s;s^>. 



(12) 



According to conformal field theory we may expect 2k p 
and 4k p CDW fluctuations (2k p = nir) in the charge sec- 
tor and 2k p SDW fluctuations in the spin sector to dom- 
inate. The onsite repulsion enhances 4kp CDW fluctua- 
tions, whereas the Heisenberg interaction term favors an- 
tiferromagnetic alignment of spins on adjacent sites and 
thus a 2kp SDW state. On the other hand, the formation 
of bound singlet pairs enhances 2k p CDW fluctuations 5 
and favors a maximum at k = 7r in the spin sector. 
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FIG. 5. Charge structure factors of the t-J-ts model at 
a = 0, 0.4 and 1 and n = 1/2, calculated by QMC simulation 
on a lattice of L = 64, fit = 32, Ari = 0.25. 



FIG. 6. Spin structure factors of the t-J-ts model at a = 
and 1 and n = 1/2, calculated by QMC simulation on a lattice 
of L = 64, fit = 32, Art = 0.25. For comparison, the result of 
a gas of free pairs S splIl (k) = ra(l — cos(fc)) is included. 



In order to investigate the effect of switching on a for 
different values J ft, we first calculate the charge struc- 
ture factor. Figure 5 shows the results at n = 1/2 for 
a = 0, 0.4 and 1 for couplings J ft = 1 and 3. Starting 
at small J ft, we expect a cusp at k = 4kp for J/t = 1 
and possibly one at k = 2fc^ also. There is a maxi- 
mum at 4k p (= 7r) for all values of a, but the cusp at 
k = 2fc^ is hardly visible at J ft = 1, since the pref- 
actor of 2k p fluctuations is much smaller than for 4k p 
fluctuations. 8 Increasing a enhances pairing of electrons 
and produces stronger 2kp fluctuations. Therefore the 
cusps at k — 2k p are much clearer for larger a. Pairing 
is also enhanced by increasing the coupling J ft. At the 
same time, a spin gap opens for the larger values of a. 
This crossover to a LEL is accompanied by a jump in the 
exponent of the 2k p CDW correlations from 1 + K p to 
K p . Accordingly there is a peak at k = 2k p diverging 
with the lattice size L in Fig. 5 for a = 1 and J ft = 3, 
while the 4kp fluctuations are strongly suppressed. 

To analyze the effect of the pair hopping term in the 
spin sector, we show in Fig. 6 the spin structure factor 
at n = 1/2 for a = and 1 for couplings J ft = 1, 3 
and 5. There is a complete spin degeneracy at J = 0, 
which is lifted by an infinitesimal value of J ft. The 
expected 2k p SDW fluctuations can be seen by a clear 
peak at k = 2kp at J ft = 1 for both values of a. With 
the enhanced formation of pairs by increasing J ft, this 
peak is completely smeared out into a broad maximum at 
k = 7r, reflecting short range nearest neighbor antiferro- 
magnetic fluctuations. For comparison, we have included 
the structure factor of a gas of free nearest neighbor pairs 
S s -pi n {k) = n(l — cos(fc)). There is a good agreement for 
a = 1 and J ft = 5 with the calculated values. 

In the t-J model (that is a = 0), the dominant an- 
tiferromagnetic configuration is a squeezed spin chain. 



The pair hopping term on the other hand favors a local 
singlet character of the spins. To investigate which ef- 
fect is dominating for a ^ 0, we calculate the correlation 
of the two spins on both sides of a single hole C s h s = 
{Sfhi + iSf +2 ) I ' (riihi + ini + 2) , and that of the two spins 
preceding a hole C ssh = {Sf S z i+l h i+2 ) / {nin i+ ihi +2 ) ■ The 
results are shown in Fig. 7 as a function of J ft for a = 
and 1 with n = 1/2 and 3/4. Since the charge and spin 
degrees of freedom are decoupled in TLL, we may ex- 
pect approximately the same values for C s h s , C ss h and 
{S\ S? +1 ). We can see in Fig. 7 that this is roughly 
fulfilled for a = 0, as C s h s and C ss h are of the same 
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FIG. 7. Three-site correlation of a hole and two spins 
Cshs = {Sfh, +1 Sf +2 )/{n,h, +1 n, + 2) and C ss h = (Sf Sf +1 h, +2 ) / 
(n l n l -yih l -\.2) in the t-J and t-J-ts model for different values 
of a, calculated by exact diagonalization of a lattice of length 
L = 16 and CSBC. 
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FIG. 8. Pair structure factor 5 pa i r (fc = 0) for the t-J-ts 
model at n = 0.5, calculated by exact diagonalization of a 
lattice of length L = 16 and CSBC. 

order of magnitude. The phase separation boundary, 
J/t ~ 3.3, for a = and n = 1/2 can be seen as a cusp 
in C s hs- On the other hand, we may expect C s hs = 
and C ss h = —1/4 for a gas of free singlet bound pairs. 
Figure 7 shows that C s hs — > and C ss h — > —1/4 with 
increasing J/t. This is due to the enhanced pairing for 
a = 1. We can see that C ss h — > — 1/4 much faster than 
Cshs — > 0. This means that next to a hole there is al- 
most always a singlet pair allowing a larger kinetic en- 
ergy by pair hopping, but due to the other terms there 
is still some antiferromagnetic correlation between two 
spins separated by a hole. 

VII. SUPERCONDUCTIVITY 

Finally we concentrate on details of the superconduct- 
ing correlations and investigate the suppression of super- 
conductivity in the low hole density region. The singlet 
superconducting structure factor is defined as 

1 L 

S pml (k) = I J2^- m HPlPm), (13) 

j,m 

where again Pj = ( c ]-|- c ] + ij, — c }i c } + i-f)/V^ is a nearest 
neighbor singlet pair. 

As pointed out in Sec. V, dominating superconduct- 
ing correlations are expected for K p > 1. We can look 
for signs of enhancement of superconductivity by calcu- 
lating the uniform component, S , pa i r (0), which contains 
both short and long range correlations. In Fig. 8 we have 
plotted S'p a i r (0) as a function of J/t for n = 1/2 and 
a = 0, 0.4 and 1. For small to intermediate couplings 
J/t < 4 the pair hopping term enhances superconduc- 
tivity, since S , pa i r (0) has larger values with increasing a. 
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FIG. 9. Real space pairing correlations for the t-J-ts model 
at a = and 1. The density is (a) n = 1/2 and (b) n = 3/4. 
Calculated by exact diagonalization of a lattice of length 
L = 16 and CSBC. 

For very strong couplings J/t > 4 an intermediate value 
of a = 0.4 is best among the three considered, producing 
the largest values of S , pa i r (0). As there is no long range 
fluctuations of superconductivity in the phase separated 
region, we can see a sharp drop in S , pa i r (0) at the bound- 
ary of phase separation, J/t ~ 3.5, for a = 0. We have 
obtained signs of enhanced superconductivity in the same 
parameter regions indirectly, through the values obtained 
for the correlation exponent K p in Sec. V. 

To separate the long range part of S , pa i r (0) we need to 
examine correlations in real space, V(r) = {P} P§) . We 
show V(r) for a = and 1 for J/t =1,3 and 5 in Fig. 9. 
For J/t = 1, superconductivity is not a dominant corre- 
lation and V(r) vanishes for the larger values of r. On 
the other hand, in the superconducting region (J/t = 3 
for a = and J/t = 3, 5 for a = 1) the correlations 
V(r) are quite large for all r. The enhancement of super- 
conductivity by the pair hopping term is demonstrated 
by the long range correlations of V(r), which are much 
larger for a = 1 than for a = 0. The coupling J/t = 5 is 
in the phase separated region for a = 0. Therefore V(r) 
is only large for r = and 8, where a pair of electrons is 
annihilated and created at either end of the Heisenberg 
chain. 

The suppression of superconductivity for high elec- 
tron densities, already discussed in Sec. V, can also be 
seen in the real space pairing correlations V(r) Fig. 9 
for n = 3/4. As for n = 1/2, the long range parts of 
V(r) are larger in the superconducting parameter region 
( J/t = 3 for a = and J/t = 10 for a = 1) than in the 
non-superconducting region. But also these long range 
parts of V(r) are much smaller compared to n = 1/2. 
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Obviously most of the contributions to S , pa i r (0) originate 
from short range correlations and superconductivity is 
suppressed. The reason for the suppression of the long 
range fluctuations {P} P§) is the lack of pairs of holes. 
The pair hopping term produces an effective repulsion 
among holes, as it has been pointed out in Sees. Ill and 
V. Therefore most of the holes are isolated, and the cor- 
relation p s h = (n;/i;+in;+2)/(l — n) in Tab. I remains 
close to one for a = 1. A strong local singlet nature of 
electrons allows the annihilation of a spin singlet pair of 
two neighboring electrons, but it is difficult to find the 
pair of holes necessary to create another nearest neigh- 
bor singlet electron pair, at a distance r ^ or 1. This is 
the main reason of suppression of superconducting cor- 
relations even though the pair hopping term enhances a 
singlet correlation between neighboring electrons. 

VIII. CONCLUSIONS 

In this paper we have numerically investigated various 
aspects of the pair hopping terms in the ID t-J-ts model 
and compared to the Hubbard and t-J models. The pair 
hopping terms derived from the Hubbard model by the 
strong coupling expansion in U jt are usually neglected 
in the t-J model. It is found that upon the introduction 
of these terms, the qualitative similarity of the Hubbard 
and t-J models is improved in the small J ft limit and 
the values of the correlation exponent K p of the t-J-ts 
and Hubbard models become similar also quantitatively 
for J/t<l. 

The pair hopping terms might be particularly impor- 
tant to superconducting correlations, since singlet elec- 
tron pairs can hop without breaking using the three site 
term, therefore a larger region of dominant superconduc- 
tivity fluctuations is expected. Indeed, our numerical 
calculations show a complete suppression of phase sep- 
aration and an extended superconducting region up to 
J ft — > oo. But the pair hopping terms also favor the 
propagation of single holes by the larger gain of kinetic 
energy compared to forming bound hole pairs, which pro- 
duces an effective repulsion among holes and there are 
practically no pairs of holes. Hence once a singlet elec- 
tron pair is annihilated, is is difficult to find the pair of 
holes needed to create it again at a some distance r ^ 
or 1 from the annihilation and superconductivity is sup- 
pressed in the realistic region for superconductivity in the 
cuprates, which is low hole density and small couplings 
J jt ~ 0.3-0.4. But we find K p > 1 already at lower val- 
ues of J ft compared to the t-J model in the low electron 
density region. 

Nevertheless we have to emphasize that due to the 
complete suppression of phase separation, the region of 
dominant superconducting fluctuations extends up to 
J ft — > oo, in contrast to many extended t-J models and 
the simple t-J model itself, where this region is restricted 
to a small precursor region of phase separation. 



In the single chain models which we examined here, the 
region of realistic parameters, J ft ~ 1/3, is not favorable 
for superconductivity and the introduction of the three- 
site t3-term does not change this conclusion. However 
in other cases, e.g. a 2D lattice or a two-chain ladder, 
the situation is quite different and the effect of the £3 
term in shifting the onsets of pairing and phase separa- 
tion boundaries may be much more significant. Actu- 
ally a recent investigation of the t-J model including a 
next nearest neighbor hopping term t' and the pair hop- 
ping terms (3 on a 4 x 4 lattice indicated enhancement 
of superconductivity. 27 We believe that it is important to 
test these observations on larger 2D lattices, where larger 
correlation lengths can be investigated, and specifically 
the realistic low hole doping region is of great interest. 
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